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c§ ■ Abstract 

We consider the restriction placed by the Bogoliubov inequality on the nature of the Bose- 
Einstein condensates (BECs) for interacting atoms in a spatial dimension D < 2 and in the 



presence of an external arbitrary potential, which may be a confining "box," a periodic, or a 
disordered potential. The atom-atom interaction gives rise to a (gauge invariance) symmetry- 



breaking term that places further restrictions on BECs in the form of a consistency proviso. The 
necessary condition for the existence of a BEC in D < 2 in all cases is macroscopic occupation 
of many single-particle momenta states with the origin a limit point (or accumulation point) of 
condensates. It is shown that the nature of BECs for noninteracting atoms in a disordered potential 
is precisely the same as that of BECs for interacting atoms in the absence of an external potential 
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I. INTRODUCTION 



Anderson disorder-induced localization (AL) describes the sudden transition of electron 
mobility from that of a conductor to that of an insulator owing to the disorder in the depths 
of the potential wells in an otherwise periodic potential [1]. The noninteracting electrons 
are in single-particles states and it may be that the introduction of the interaction between 
electrons may lead to derealization. This has prompted the study of AL of ultracold atoms 
in ID disordered optical potentials. The ultracold atoms are considered to be noninteracting 

nn 

and in a BEC [2-5]. This bring to the fore the fundamental question of the existence or 
absence of (gauge invariance) symmetry breaking terms in the Hamiltonian of interacting 
Bose gases in spatial dimensionality D < 2 that give rise to macroscopic occupation at finite 
temperatures of many single-particle momentum states. 

Existing proofs of the absence of a BEC for D < 2 are based on Bogoliubov's inequality 



6]. The proofs either (1) introduce a symmetry-breaking fie 
to get a non-vanishing order parameter at finite volume 
the static order-order correlation function (with or without a symmetry-breaking field) [12} . 



d into the Hamiltonian in order 
or (2) work directly with 



However, none of these proofs is completely rigorous [13] since they do not take into account 
the infinite-dimensional nature of the Hilbert space of a Bose system in a finite volume and 
the consequent unbounded nature of the operators that appear in the Bogoliubov inequality. 
Thus, the validity of the Bogoliubov inequality itself is not established beforehand, which 
includes the proper h andlmg of the th^od^c « ,3- -far a S the results of 
the proofs are concerned either (1) a BEC in the single-particle state with momentum p, 
usually p = 0, is proven to be absent, or (2) a BEC in any single-particle state is shown 



to be absent. The former result has been established 13] completely rigorously; while in 
deriving the latter results jl-ioj], the Bogoliubov inequality is used without first proving it 
specifically for the infinite-dimensional case of the Bose gas. 

II. INTERACTING BOSE GAS 

Consider the Hamiltonian for an interacting Bose gas 

H = J rfr^ t (r)(^-V 2 )^(r) + J drft '(r)V ext (r)ip(r) + J drdr'ft(r)ft( r >)V(r-r')4>(r')4)(r), 

(1) 
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where V ex t(r) is the external potential, V(r — r') is the two-particle, local interaction po- 
tential, and V'(r) and ^(r) are bosonic field operators that destroy or create a particle 
at spatial position r. For the case of periodic potentials, we do not consider interactions 
between the bosons and the crystal that result in the creation of phonons. Macroscopic oc- 
cupation in the single-particle state ip(r) result in the no n- vanishing [ijj] of the quasi- average 
ip(r) =< V>(r) > and so the boson field operator 

^(r)=V(r) + 0(r), (2) 

with 
and 

El£ k 'l 2 = i, (4) 

k' 

where N is the number of atoms in the condensate and V(D) is the D-dimensional "volume" 
and < 0(r) >= 0. The operator <p(r) has no single-particle states that are in the condensate 
and so / drip' (r)-0(r) = 0. The separation of ^(r) into two parts gives rise to the following 
(gauge invariance) symmetry breaking term in the Hamiltonian (1) 

Hsymm = J dv^ (v)^(v) J dv'[V{v - v') + V{v' - v)]\^{v')\ 2 + h.C. 

= J dr(p\r)x(r) + h.c. (5) 

The presence of this nonzero H symm in the Hamiltonian gives rise to further macroscopic 
occupation in states other than the original state given by ip(v) and so the condensate wave- 
function ^(r) gets modified by augmenting the single-particles states where macroscopic 
occupation occurs. In such a case, macroscopic occupation in the state b would give rise 
to macroscopic occupation in the states a, such that a ^ b, whenever the matrix element 
< ab\V\bb > of the potential V, which is the last term in Eq. (1), does not vanish. For 
noninteracting particles in an external potential, the Hamiltonian is diagonal in the repre- 
sentation of the eigenstates of the Hamiltonian and so macroscopic occupation in any given 
eigenstate of the Hamiltonian does not give rise to further macroscopic occupation in any 
other energy eigenstate. However, the effect of interparticle interactions can generate macro- 
scopic occupation in other energy eigenstates. For instance, for a harmonic trap, macroscopic 
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occupation of the lowest energy state, the ground state, does not give rise to macroscopic 
occupation of any of the higher energy eigenstates. However, interparticle interactions may 
allow the generation of macroscopic occupation in other energy eigenstate thus modifying 
the original BEC. Note, however, that macroscopic occupation only in the single-particle 
state with momentum p does not give rise to macroscopic occupation in any other momen- 
tum state since the matrix element in the momentum representation < qp|V|pp > vanishes 
by momentum conservation unless q = p. 

This consistency proviso requires that the correct condensate wavefunction if>(r) corre- 
sponds to that which gives rise to no symmetry breaking term in the Hamiltonian. That is 
to say, Hgymm vanishes for the correct condensate wavefunction ip(r). For instance, macro- 



scopic occupation in the single-particle states with momenta 0, ±k gives rise [15j, with the 
aid of the symmetry breaking term H symm and owing to linear momentum conservation, 
to macroscopic occupation in the single-particle momenta states ±2k ; therefore, </?'(r) is 
orthogonal to both ijj(r) and x(r) and thus one has the possibility of macroscopic occupation 
in all the momentum states 0, ±ko, ±2ko, ±3ko, • • ■■ In general, one can have also a finite 
or an infinite sum over the momentum variable ko and so the momentum conserving 
two-particle interaction in the Hamiltonian (1) allows for condensates in the single-particle 
momentum states k + J2i m ikj, where rrii = 0, ±1, ±2, • • • [16[. The momenta {kj} represent 
a set of vectors that are, in general, incommensurate thus giving rise to nonperiodic or ape- 
riodic condensates. However, if one or more of the momenta kj approaches zero, then the 
sequence of condensates over the momenta sequence {k'} with nonzero £i> has k' = as a 
limit point (or accumulation point) as kj — > thus removing the l/fc 2 -singularity 17| that is 
necessary in all proofs of the absence of a BEC for D < 2. For a system confined to a box of 
length L by an external potential, momenta is quantized and so k = ^n x x+ ^-^y + ^n z z. 

Note, however, that recent papers still suppose that interacting Bose gases in D < 2 and 
at T > possess no BECs 18M20 . The phase transition is described 18[ by the Berezinskii- 
Kosterlitz-Thouless theory 2lM23|. which does not involve any spontaneous symmetry break- 
ing but, instead, is associated with a topological order embodied in the pairing of vortices 
with opposite circulations. It is interesting that both the singly quantized vortex in a doubly 
connected system and the condensate that gives rise to two vortices of equal and opposite 
circulations are both described by condensates with a macroscopic occupation of infinitely 
many single-particle states with fc = 0a limit-point Similarly, limits on Bose-Einstein 
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condensation in confined solid 4 He, where large superfluid fractions have been reported, is 



24J on supposing macroscopic occupation in only the single-particle momentum state 

C. In addition, it is supposed that for 



based 

k = 0, which is not the most general, possible B 



2D, superfluidity is not a consequence of a BEC |21l . |22| but is associated with the onset 
of algebraically decaying off-diagonal long-range order (ODLRO) 25J. However, the exis- 
tence of ODLRO can be a consequence of macroscopic occupation in many single-particle 
momenta states, which is equivalent to the existence of a BEC. In what follows, we con- 
sider macroscopic occupation for D < 2 in many-momenta states with k = a point of 
accumulation. 



A. Example of BEC in ID 

Consider the following one- dimensional example for a BEC with macroscopic occupation 
in the single-particle momentum states k = k^n with n = 0, ±1, ±2, • • • for a Bose system in 
the presence of a nonperiodic potential, where k = 2tt/L. The BEC (3) becomes 

,/ x [No ( k x-i/2 +2? k cos{k nx) 

where v > so that ip(0) is finite. The BEC is normalized as follows J 2w ^ k ° dx\ijj(x)\ 2 = Nq. 
Expansion (6) represents a Fourier series and so tp(x + 2n /ko) = ip(x), which is a dynamical 
consequence since macroscopic occupation in the single-particle momentum states k = 0, ±/c 
implies macroscopic occupation in the states with momenta k = nko with n = 0, ±1, ±2, • • • 
owing to the symmetry-breaking term (5). As a function of the complex variable ko, the 
Fourier series in (6) possesses simple or higher order poles at ko = ±itz/n, n = 1,2,3, • • •, 
for v + 1/2 = positive integer or branch point singularities at ko = ±in/n, n — 1, 2, 3, • - -, 
for u+ 1/2 7^ positive integer. Therefore, ko = is an accumulation point of poles or branch 
point singularities depending on the value of v. 

The sum (6) can be approximated with great accuracy by an integral when ko is arbitrarily 
small. This would represent the passage of the Fourier series for the periodic function to 
a Riemann integral for a nonperiodic function. Accordingly, in the limit k — > 0, k = 
becomes a point of accumulation of condensates, the sum (6) may be converted into an 
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integral, and so 




(k 2 + K 2 ) 2v + l 



1 



-1/2 f+°° 



— oo 



dA; 



(k 2 + « 2 )-+V2 



cos{kx) 



1 «r(2i/ + l)iV 



k v \x\ u K u (k\x\), (— oo < rr < oo), 



(7) 



2^-1/2^1/4^ r(2z/ + l/2)r 2 (z/ + l/2) 



where T(z) is the gamma function and K v {z) is the modified Bessel function of the second 
kind. The integral (7) converges for \x\ > provided v > —1/2; however, if the condensate 
wavefunction is required to be bounded at x — 0, then (7) converges for \x\ > provided 
v > 0. 

Now d[^^(^)]/dz = -z u K u ^ 1 (z), K v (z) ->■ 2 1 ^ 1 r(z/)/^ as * ->■ for fti/ > 0, and 
= K v (z). Therefore, z u K v (z) = 2 v ' 1 T(u) - 2 iy " 3 r(z/ - l)z 2 + • • • as z ->■ for 
z/ > 1 and so we have a local maximum at the middle of the BEC given by (7). Next one 
has that z"K v (z) = 2 !y ~ 1 r(z/) - T(l - v)z 2v / \2 u+1 v) + • • • as z ->■ for < v < 1 and so 
for < v < 1/2 the BEC (7) has a cusp singularity at re = with the derivative tending 
toward oo or — oo as one approaches the cusp. For v = 1/2, the BEC (7) is continuous 



at x = but the derivative is discontinuous there with slope of ^/tt/2 or -y/vr/2 as one 
approaches x — 0. For 1/2 < v < 1, the BEC (7) attains its largest value of the center 
of the localized BEC but it is not a local maximum. The extremal cases of v — 0, 1 are 
as follows. For v — 0, -^(x) oc K (k\x\) = — In |x| + ••• thus there is a cusp singularity 
with slope oo or — oo as one traverses x — 0. Notice that BEC (7) is not bounded at the 
center of the localized condensate for v — 0; nonetheless, the integral of the BEC density, 
which gives the total number of particles in the condensate, is finite. Finally, for v — 1, 
zK 1 (z) = 1 + (z 2 /2) \nz + • • • as z — > and so -^(rr) attains its largest value at z = but it 
is not a local maximum. 

The standard deviations Ax and Ap follow directly from (7) 







and so 



(Ax)(Ap) 



h (u+ l/4)(z/+ 1/2) 1 
2^ (*/ + !)(*/ -1/4) " 2 



(9) 



for ^ > 1/4. 
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The sum in (6) can be carried out explicitly for v — 1/2 and one obtains 



where (f>(x + 2n/k ) = <f>(x) and the normalization constant A is 





(10) 



2 



(11) 



Note that for ko — > 0, one has that 



ip(x) = \J kNq e (-co < rr < oo), 



(12) 



which agrees with Eq. (7) since K 1/2 (z) = ^/2z e~ z . Result (12) corresponds to a localized 
BEC that decays exponentially and has a discontinuous derivative at the center. 

For v > 0, K u (z) — > 2 u ~ 1 T(v) / ' z v as z — > and so the condensate wavefunction ip(x) 



given by (7) is finite at x — 0. In addition, K v {z) —¥ \j^/2z e~ z as z — >■ oo for z/ > —1/2 

and so the condensate wavefunction (7) is similarly localized. Note that prior to taking 
the limit ko — > 0, which approximates the sum (6) with the integral (7), the condensate 
wavefunction is periodic, viz., ip(x + 2n/ko) = ip(x). However, in the limit ko — > 0, one has 
an even, nonperiodic, localized condensate wavefunction. 

Actually, one can interchange the roles of x and k in Eq. (7) and so one has, instead, 
a condensate wavefunction ip(x) oc (x 2 + a 2 )^ u ' 1 ^ 2 with a corresponding momentum distri- 
bution (p(k) oc Ik]" K u (a\k\) , where a is a length scale. It is interesting that in experiments 
with ultracold atoms (Ref. 28), BECs have been found that suggest a power-law decrease in 
the wings of the atomic density with an exponent close to the value 2, viz., |-?/>(a;)| 2 oc x~ 2 as 
\x\ — > oo. This would correspond to v — and so | 2 oc (x 2 + a 2 )' 1 for (-co < x < oo). 
The momentum distribution associated with this BEC density is given by the modified 
Bessel function of order zero (p(k) oc K (a\k\) for (— oo < k < oo), where the singular point 
k = is a logarithmic branch point since K (z) oc — ln^ as z — > 0. The standard deviation 
Ap = h/(2aV2) is finite even though <p(k) diverges logarithmically as k — > 0; however, Ax 
is infinite even though ^(x)) 2 is bounded everywhere. Note that the behavior of ip{x) for 
x ^> a is determined by the behavior of its Fourier transform (p(k) for ka 1. 

It is interesting to consider a finite rather than an infinite sum in Eq. (6), 





( £ 



n=-N/2 



N/2 



(kin 2 + k 2 ) 2 ^ 1 




=-N/2 



(kin 2 + K 2 )-+i/2' 



k cos(k nx) 



(0<x< 2n/k ), 



(13) 
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and thus consider the combined limit k — > and N — > oo in order to approximate the sum 
by an integral. If k$N — > oo as k — > and N — >• oo, then one obtains the previous result 
given by Eq. (7). However, if k N/2 — > K < oo as k — > and N — > oo, then 

. /iVo / r K 1 x-i/2 cos(xA;) 

Y I = 1 — / d/c--T „. „ — / Cik-—: — — 7T7o i — OO < X < oo . 14 

AW V 2tt W-x (k 2 + k, 2 ) 2u+1 J J~k (k 2 + K 2 Y +l / 2 ' y i \ ) 

For one-dimensional crystals, K < n/a, where a is the length of a unit cell. The condensate 
wavefunction x( x ) is normalized by f^dx\x(x)\ 2 = N . Note that for K\x\ 3> 1, the 
second integral in (14) is dominated by the small values of k, viz., k K, and so the 
range of the integral in (14) can be extended to ±oo when K\x\ ^> 1 and one has that 
x(x) oc ^^^(/tlxl) for K\x\ ^> 1, which becomes x( x ) °t l^l^ -1 ^ 2 e~ K ' :r ' for k\x\ ^> 1. 



B. BEC in the harmonic trap 

If the external potential is a one-dimensional harmonic oscillator, then the field operator 
ip(x) is expanded in terms of the energy eigenstates ip n {x) of the harmonic oscillator. If one 
has macroscopic occupation in the ground state ^po(x), then ^p(x) = ipo{x) + <p(x), where 
the operator (p{x) is orthogonal to ipo(x) and so <p(x) has nonzero expansion coefficients for 
only the creation operators of the higher energy harmonic oscillator eigenstates. It should 
be noted that if one considers atom-atom interactions, then it may be that states other than 
the ground state may be also macroscopically occupied. The latter occurs if the matrix 
element < ZO | 1 00 > of the atom- atom interaction potential V does not vanish for I ^ in 
which case there would be additional macroscopic occupations in the higher energy harmonic 
oscillator eigenstates ipi(x). 

Consider the following illustration of a condensate in the ground state of a one- 
dimensional harmonic trap 

M*) = \hf ( E kve-^y 1 ' 2 f; A; e-^ 2 cos(A; na:), -n/k < x < ir/k . 

' n=— oo n=— oo 

(15) 

The condensate (15) has macroscopic occupation in the single-particle momentum states 
kon, with n = 0, ±1,±2, ip(x) is periodic ip(x + 2n/ko) = ip(x), and is normalized 
as follows Jll k ° ko dx\ijj(x)\ 2 = N . Each Fourier coefficient of the series (15) possesses an 
essential singularity at k = oo and in the limit ko — > 0, may be approximated arbitrarily 
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well by a Riemann integral and so 

Mx) = M( r dke-^ k2 Y 1/2 r dk e^ k2 cos(kx) 

V 2lT ^ J-oo ' J-oo 

= ^(^) V V^, -oo<x<oo. (16) 

Therefore, the condensate associated with the lowest energy state in the harmonic trap 
is represented by the macroscopic occupation of single-particle momentum states given by 
k = nko, where n — 0, ±1, ±2, • • -, with fc = 0a point of accumulation as ko — > 0. 

Note that as a function of the complex variable ko, the BEC given by the series (15) con- 
verges for K/cq = (Jftk ) 2 — (Sfco) 2 > and diverges for K/cq < 0. The regions of convergence 
and divergence resemble a Minkowski spacetime diagram. The series (15) converges in the 
spacelike regions and diverges in the timelike regions and the light cones. Along the real 
axis, the series (15) converges except the singularity that it encounters at the origin k = 0. 



III. BOGOLIUBOV INEQUALITY 

The absence or presence of a BEC in spatial dimensions D < 2 is based on Bogoliubov's 
inequality 

\({A,ti}) >k B T\([C,A])\*/([[C,H],Cl]), (17) 

where H is the Hamiltonian (1) of the system with arbitrary local interparticle and external 
potentials, the brackets denote thermal averages, and the operators A and C are arbitrary 
provided all averages exist. 



Consider the following operators 



13, 



C = J c/re lk - r ^ t (r)^(r) (18) 



and 



A = Jdrj dv'e- lk - r f(r)f*(r')4>\r)4>(r'), (19) 
where k is arbitrary. Now, 

([At,A]) = V(D)([(7,i]) = y 2 (D){A -A |A k | 2 -^(a^ q )|e q+k | 2 }, (20) 

q 

([[C,H],&)) = ^f N , (21) 

fit 



with 

^k = E^ +k = ^/dr|^(r)|V k - r (22) 

with the aid of Eq. (3), where N is the total number of particles. Now, \A-^\ < 1 by the 
Cauchy-Schwarz inequality where the equality holds for k = 0, that is, A — 1 with the 
aid of Eq. (4). The vector q ^ {k'}, where {k'} is the set of condensate vectors for which 
f k , jL 0. For a BEC at rest, |£ k | 2 = |£_ k | 2 . Note that £ q+k ^ for (q + k) G {k'} and 
£q+k = for q ^ {k'} and k G {k'}. 

We sum the Bogoliubov inequality (17) over the single-particle momentum states in the 
set {k'} constituting the condensate, which includes an arbitrary neighborhood of the point 
of accumulation of the condensate at k' = that corresponds to a condensate at rest. We 
want to find an upper bound of the anticommutator ({A, A*}) = 2(A4 t ) + ([A^,A]). We 
extend the sum over the first term (A<4.t) over all values of k thus obtaining a larger upper 
bound 

MN NV 2 (D) > £(AA f >, (23) 
k 

where the condensate wavefunction ipij) is orthogonal to the operator <£t(r) and we assume 
that the condensate density is bounded from above by |/(r)| 2 < M with 1 < M < oo 
since / <ir|/(r)| 2 = V(D). In (23) use has been made of the completeness relation for the 
momentum eigenstates and a negative term resulting from a single commutation has been 
dropped. Note that we are considering a condensate where all the single-particle states 
with momentum k' are occupied macroscopically with k' = a point of accumulation. In 
addition, we are supposing that the number of particles in the "volume" V(D) is fixed, that 
is, we are employing a canonical ensemble and so / drfo (r)i[j(r) = J] k ^ k ^ k — N is actually 
the c-number N. 

Consider next the sum over k' of the commutator ([A 1 ", A]), 

£<[iU]> = W 2 P)£[iH^ k f], (24) 
k' k' 

with the aid of (20) and where £q +k = for q ^ {k'} and k G {k'}. This sum over the 
commutator is bounded from above provided the sum is restricted to values of k' that have 
a finite, upper bound. Now the right-hand side (RHS) of inequality (17) becomes 

k' 
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with the aid of Eqs. (20) and (21). Note that the RHS is bounded in the upper limit of the 
sum; however, it is the lower limit as k' — > for D < 2 where the sum may diverge which 
would result in no BECs, viz., N = for T > 0. Combining Eqs. (23)-(25), we have for 
the Bogoliubov inequality, 

M^^O-W)^^^ (26, 
k k 

Therefore, the existence of a BEC for T > requires the convergence of the sum on the 

RHS of (26) over the macroscopically occupied single-particle momentum states k' of the 

condensate. Note that the sums in (26) over the condensate momenta can be approximated 

by integrals according to J2\r> — > V(D) J dk' and so 

The integral on the RHS has no infrared divergence since by (22), (1 — l^^'l 2 ) vanishes 
quadratically as k' — > and so the l//c 2 -singularity is removed thus allowing the existence 
of a BEC for D < 2. 



A. Removal of 1/k 2 singularity in ID 

To illustrate the removal of the 1 //^-singularity for the existence of a BEC for D < 2, 
consider the ID example given by Eq. (6) where 

= (kW (» = 0,±1,±2,-), (28) 

with the normalization constant B given by 



S " ( f £o(*8n 2 + « 2 ^ l) " 1/2 ^ 



with the aid of (4) and so by (22) 

2 +oc 1 1 

A kol = A. kol = B -^————j- _______ __. (30) 

Now 

r ^ +1 / 2 ) - r dx e -[*g(n-0 a +« a ]« 

[jfe2( n _Z)2 + «2] V+ l/2 J 
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" . r(j + ^+l/2)^(/ 2 -2n/)^ 

j ] -{kW + «2)i+H-l/2 lJ J 

for ^ > —1/2, which gives the following series expansion for (30) in the neighborhood of 
Z = 0, 

2 I> + 3/2) +2? fc 2 Z 2 
r(z, + l/2) n £- 0O (^ + ^)^ 2 + 

I 0.2 ^ + 5/2) +2? fcgZV 4 
The sum on the RHS of (26) is 

+oo 1 

E p»(l-4w) 2 (l + 4w) a (33) 

and so the singularity at / = in the summand is removed as indicated by the expansion 
(32) of A kol . 

Similarly for the BEC in the harmonic trap where (15) gives that 

oo 1 oo 

4w=( E e- 2 ^ 2 " 2 )- 1 e-^ 2 e-^ l \oM^lkl), (34) 



n=— oo 



which on expanding in powers of I about I = removes the /^^-singularity at k l = in 
the series on the RHS of the Bogoliubov inequality (26). 

IV. BECS IN PERIODIC AND DISORDERED POTENTIALS 

The existence of su per fluidity is usually associated with the existence of a BEC. The 



existence of superflow 26| in solid helium 4 He has stimulated the search of a BEC in solid 
helium thus establishing the existence of BECs in all three states of matter-gas, liquid, and 
solid. It is to be noted that the proofs of the absence of a BEC in gases and liquids for 
D < 2 applies also for atoms in external infinite periodic potentials and thus to crystalline 
solids. If one supposes macroscopic occupation in the momenta states k, k ± q , then the 
symmetry breaking term (5) gives rise to macroscopic occupation in the momenta states 
k ± nq with n = 0, ±1, ±2, • • •. Accordingly, the condensate wavefunction is given by 



= V W) Sj^"^"^ S e ' k "' ^ ' k(r) • <35> 
which is of the Bloch form since W] £ (r) has the periodicity of the lattice, that is, u-^(r) = 
u^(r+t m ) since e l( ^o"^™ = 1. The primitive lattice translation vector t m = mia+m2b + m3C, 
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where rrii can take all integer values and a, b, and c are the edges of the unit cell, which forms 
a parallelepiped. Therefore, the possible values of the lattice condensate vector(s) q are 
in q-space or reciprocal space. Note that Bloch functions ^(r) and ^'(r) are orthogonal 
for k' 7^ k provided |k| < |q |/2 and |k'| < |q |/2 and so / dr^j^r)^ (r) = N 8^' ^ . The 
expectation value of the momentum in the condensate follows from (35) and so 

/oo 
dr^(r)(-^V)V k (r) = N h £ (k + nq )|£ k+nq J 2 = N hk (36) 
n=— oo 

provided |Ck+ n q 1 2 = l^k-nq I 2, Therefore, k represents the energy dependent Bloch wave 
number or the quasi momentum per particle Kk of the condensate. 

For the case of a ID crystal of length L = Na, where there are N primitive cells of length 
a one has that 

fa(x) = e ikx u k (x), (37) 

with Uk(x + a) = Uk(x), where k = 2rrm/Na (m = 0, ±1, ±2, • • •) owing to the boundary 
condition ip k (x + N a) = ipk( x )- The range for the momentum k is given by —n/a < k < n/a 
since if the Bloch condition holds for k it also holds for k' = k + 27im/a. One has from the 
Schrodinger equation that ^p-k(x) = ipl(x) and so the negative values of k do not give rise 
to new solutions. In addition, the values of m are limited to m = 0,1,2, N — 1 since 
the higher values of m do not generate any new solutions and so we have only iV solutions 
per band. In fact, one needs an additional label, a band index, to be added to the Bloch 
function to describe which band the function belongs. 

A. Atom-atom interactions and disordered potentials 

The Bloch form for the condensate wavefunction ^(r) given by (35) is appropriate for 
noninteracting particles in a perfect crystal for D = 3. However, the condensate (35) does 
not remove the 1 //^-singularity for finite q for spatial dimensions D < 2 and T > 0, which 
is required for the existence of a BEC for systems of noninteracting or interacting particles 
when embedded in periodic or in disordered potentials. It is interesting that the removal 
of the 1 //^-singularity leads to localization. Therefore, we suppose that the inclusion of an 
interparticle potential, given by the third term in the RHS of (1), gives rise to a BEC provided 
q is arbitrarily small and so (35) ceases to be of the Bloch form. It is interesting that this 
behavior is equivalent to supposing a BEC that is a linear superposition of condensates with 
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differing discrete, translational motion and so the resulting condensate wavefunction is given 
by 

VW = EV*" r «kOO> (38) 
k 

where the prime in the sum indicates that |k| < | q 1 / 2 . The existence of a point of accumu- 
lation, or limit-point, at k = is what is required for the existence of a BEC for systems with 
spatial dimensionality D < 2. Expression (38) represents a sort of wave packet constructed 
not by means of plane waves but instead by means of plane waves modulated by a Bloch 
state, which still leads to an expansion in terms of plane waves as the Fourier series given by 
(3). It is interesting that such types of states have been used in the study of the dynamics 



of electron wave packets in crystals 27|. In addition, the wave packets so prepared remain 



in the same band at later times, in our case the lowest energy band, and are referred to as 
Bloch-type states [27] since even though they are not Bloch states, viz. there is no single 
momentum ftk such that ip(v + t m ) = e*^'^ m ^(r) even though « k (r + t m ) = u k (r). The 
momentum associated with the condensate (38) is 

f drf(r)(- 2 W)^(r) = iVo^'k|a k | 2 = 0, (39) 

k 

where the average momentum of the condensate is zero, that is, we choose the system with 
respect to which our condensate of N particles is at rest and so |a k | 2 = |ct_ k | 2 . Similarly, 
the kinetic energy T associated with the condensate (38) is 



(40) 



2/i /rv ' 2/i 1 2/i k 



where /1 is the particle mass and use has been made of the normalization conditions 

E k > k | 2 = land E^_ooie k+ „ qo l 2 = l- 

One would expect that the kinetic energy T of the condensate should be rather small, 
which suggest neglecting the second term on the RHS of (40). This corresponds to replacing 
w k (r) in (38) by a constant. Accordingly, 

^(r)^£'«k elkT (41) 
k 

Note that the BEC (41) is not a Bloch function since ^(r + t m ) 7^ "^(r). In addition, for 
cases where the sum over k is associated with a nonisolated singularity at k = 0, such a BEC 
is equally applicable to Anderson localization of ultracold atoms in a disordered potential. 
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B. BECs in ID crystals 



Consider the BEC in the lowest, or first Brillouin zone, energy band with macro- 
scopic occupation in the first, say, 2M + 1, lowest energy states, viz. with k = 
0, ±27c/Na, ±Arr/Na, ±2irM/Na with < M < N/2. Thus, the BEC (41) becomes 

M 

1>(x) = £ » m e"/ Na , (42) 

m=~M 

where a m = a_ m with normalization condition ^ = J2m=-M \ a m\ 2 - Therefore, in the 
domain \ka\ < 2irM/N, our condensate possesses macroscopic occupation of 2M single- 
particle momentum states. Now ip(x) in (42) satisfies the periodic boundary condition 
■ip(x + L) = ip(x); however, ip(x) does not satisfy the Bloch condition, viz., ip(x + la) ^ ip(x) 
for < / < N. Notice that in the limit N — >• oo, M — >• oo such that M/N -> K, one has, in 
any arbitrary neighborhood of the lowest energy state with k = 0, an unlimited number of 
nonvanishing a m with k = an nonisolated singularity. This behavior is what is required 
to remove the l/A; 2 -singularity that appears in the Bogoliubov inequality in order for the 
crystalline system to possess a BEC in ID. 

C. BECs in disordered potentials 

The general expression of BECs in a periodic potential is given by 
*k( r » = S £ «k + ,„q 1+w q 2+ ... e i|k+ '"1- + "* + >W kr 1% (r), (43) 

V V / ni,ri2,---=— oo 

with ^^(r) = u-^(r + t m ) for any primitive lattice translation vector t m and where the 
set of vectors {q { } are in the reciprocal lattice space and so e^'* m = 1 for i — 0, 1, 2, • • •. 
The generation of macroscopic occupation in the momenta states in Eq. (43), which occurs, 
albeit, even in the presence of arbitrarily weak two-body interactions, is a direct consequence 
of the symmetry breaking term (5) and the supposition that the single-particle momenta 
states 0, ±q 1 , ±q 2 , ±q 3 , • • • are macroscopically occupied. Of course, if the interparticle 
potentials are not negligible, then the BEC cannot be of the Bloch form (43). Therefore, 
in the presence of interparticle interactions, not all the vectors {q^} are in the reciprocal 
lattice space, in which case the BEC can still be expressed in the form of Eq. (43) except 
that now ^^(r) ^ + t m ) and so the BEC is not of the Bloch form. The latter is what 
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one would expect for nonperiodic or disordered potentials for spatial dimensions D = 3. It 
is clear that for D < 2, some of the vectors {q^} must vanish in a limiting process in order 
to remove the l//c 2 -singularity in the Bogoliubov inequality (26) required for the existence 
of a BEC in spatial dimension D < 2. 

It is important to remark that the existence of a BEC for D < 2, which requires the 
removal of the l//c 2 -singularity by means of a point of accumulation at k = 0, is precisely 
the same for both noninteracting atoms in a disordered potentials as well as for interacting 
atoms in the absence of an external potential. Therefore, the mathematical form of the BEC 
in D < 2 and T > for noninteracting particles in a disorder potential is indistinguishable 
from that of interacting atoms in a uniform medium. 



V. ANDERSON LOCALIZATION OF ULTRACOLD ATOMS 



Anderson localization of matter waves has been observed with cold atoms from a nonin- 
teracting BEC in a one-dimensional disordered potential generated by a laser speckle pattern 



28( and where the quasi-periodic lattice is the result of the addition of noncommensurate 



optical periods 29]. The experiment consists in releasing a BEC in the ID disordered op- 
tical potential, where all the noninteracting atoms are originally in the same single-atom 
wavefunction, viz. the localized condensate. The atomic wavefunction initially expands and 
subsequently stops expanding and the resulting wavepacket has wings that decay exponen- 



tially or as a power-law [28| and exponentially or Gaussian-like 29]. 

The condensate of a nonideal Bose gas in ID and in the absence of an external potential is 
not given by the macroscopic occupation of a single momentum state, as is the case in an ideal 
Bose gas in D = 3, owing to the l//c 2 -singularity in the Bogoliubov inequality. Actually, 
for spatial dimensions D < 2, the condensate must to be nonuniform 17( and localized as 
shown above. Therefore, both noninteracting bosons in disordered potentials and interacting 
bosons in the absence of external potentials give rise to Anderson localization. It would be 
interesting if one could experimentally discern the relative contribution to localization by 
the disordered potential and by the interactions among atoms. In particular, if localization 
would persist in an ultracold, non-dilute atomic ID gas on expansion, where atom-atom 
interactions cannot be neglected, even in the absence of a disordered potential, viz., in the 
absence of any external potential. The experimental proof that a theory, based on the 
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requirements imposed by the Bogoliubov inequality, allows for localization with repulsive 
atom-atom interactions present would lend some support to the original work of Anderson 
of a sudden phase transition from conductor to insulator via the degree of disorder in the 

n 

material [1]. It is interesting that the Mott metal-insulator transition leads to localization 
without randomness owing to electron-electron interactions, which is somewhat similar to 
localization in an interacting Bose gas for D < 2 in the absence of a disordered potential. 



Localization has been studied for expanding BECs in weak random potentials [30] . The 
localization of a single particle is treated in the Born approximation and the correspond- 
ing Lyapunov exponent, characterizing the spatial asymptotic decay of the BEC density, is 
determined by values of the correlation length of the disorder <jr and the high-momentum 
cutoff at the inverse healing length l/£j„. For £ in > a R , the BEC wave function is expo- 



nentially localized, whereas for £ in < or, the spatial decay is algebraic |30j. For the speckle 
potential considered, the Fourier transform of the correlation function vanishes for momenta 
k > 2(jft 1 resulting in a vanishing Lyapunov exponent for k > cr^ 1 in the Born approxima- 
tion. In one-dimensional elastic scattering, the particle wave vector k in forward scattering 
remains unchanged or changes sign in backscattering resulting in a momentum transfer of 
2k. Therefore, the study of localization for k > a^ 1 requires going beyond the Born ap- 
proximation. The higher order corrections to the Born approximation give rise to "effective 
mobility edges" at k = pd^ 1 , where p is an integer that characterizes the successive correc- 



tions to the Born approximation 3l|, |32] . It is interesting that these higher-order terms of 



the Born series are necessary even for k < . 

The atom-atom interaction in the Hamiltonian (1) requires macroscopic occupation in 
momenta k = nk , where n = 0,±1,±2, owing to linear momentum conservation, if 
there is macroscopic occupation in the two momenta states k = 0, ko. Therefore, for ID, 
forward scattering of momentum k requires backward scattering into the momentum state 
—k. Our example of a BEC in ID given by (6) reflects such requirements, which is a signature 
that one is dealing with a nonperturbative feature of the interparticle potential. Notice that 
in obtaining the BEC density from the BEC wavefunction (6), one does not suppose that 
the phases for different momenta are uncorrelated since such supposition would give rise to a 
constant BEC density, which would be equivalent to supposing macroscopic occupation in a 
single momentum state and thus to a uniform BEC. The assumption that localized function 
for a given momentum k are uncorrelated is made in the case of random potentials [0. 
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The sum (3) representing the condensate wave function can be approximated arbitrarily 
well by an integral in any spatial dimension. This is especially important for D < 2 where 
the BEC cannot occur in a single-particle momentum state but must be accompanied by 
a point of accumulation of single-particle momenta condensates at k = and so the BEC 
density is spatially nonuniform. For instance, the example for ID given by the sum (6) 
with its corresponding approximate value given by integral (7). For the case v = 1/2, the 
integrand of the Fourier cosine transform in (7) possesses two simple poles at k — ±m 
and so for x > (x < 0), one calculates the integral by closing the contour on the upper 
(lower) half-plane thus enclosing the simple pole at k = in [k = —in ) that yields result 
(12). Therefore, the long-tail behavior of the BEC density is determined by the singularity 
in the complex fc-space closest to the real axis. Note that for v > —1/2, the nature of 
the singularities of the integrand in (7) is generally branch points or poles of higher orders 
at k = ±m. Nonetheless, the asymptotic formula is still given by an exponential decay, 
which is a property of the Bessel function K v (z) that tends exponentially to zero as z — > oo 
through positive values. An exponentially decaying BEC density is a direct consequence of 
Fourier transform <f(k) of ip{ x ) given by a meromorphic function of k, viz., ip(k) is analytic 
except at a set of isolated points, e.g., ratios of rational functions of k. If a pole of (p(k) 
is off the imaginary axis, then the asymptotic exponential decay of ip(x) is modulated by 
sinusoidal functions, which would represent a remnant periodicity in the system. 

In Section II. A, we mentioned examples of BEC in ID with algebraic localization, viz., 

1 1 r°° 

X (** + a')** V2 = (^T> + 1/2) loo dk W K Mk\)cos(kx) (44) 

for v > —1/2, a > 0, and — oo < x < oo. The function z v K v {z) is an analytic function of 
z for v = n + 1/2, (n = 0,1,2,---), but the origin z = is a branch point singularity for 
v 7^ n + 1/2, (n = 0, 1, 2, • • ■) [33]. However, the function that appears in the integrand in 
(44), viz., l^l^i^al^D, is not an analytic function of z since \z\ is not an analytic function of 
z. In fact, K v (z) is an analytic function of z throughout the z-plane cut along the negative 
real axis since z = is a branch point singularity. According to the Riemann-Lebesgue 
lemma, the Fourier representation of the BEC ip{x) goes to zero as x — > oo and so the 
large distance behavior of the BEC is determined by the small k behavior of its Fourier 
transform ip(k). Our example (44) of algebraic decay suggests that the behavior of the 
BEC ip(x) oc l/(x 2 + a 2 ) u+1 / 2 is determined as x — > oo by a branch point singularity at 
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k = 0. Note also that the value of the BEC ip{x) for x ~ is determined by the behavior 
of its Fourier transform ^p(k) over a finite range of values of k near k = for v > since 
z v K v (z) — > 2 I/ ~ 1 T(z/) as 2 — > for z/ > 0. However, for v = 0, the range of values that 
contributes to the integral is very small and quite close to k = owing to K (z) — > — In 2 
as z — > 0. For instance, the contribution to ip(0) by the integral (44) for v = 1/2, is 63% 
from the region < & < 1/a and 37% from the region k > 1/a. On the other hand, for 
v = 0, the contribution to ip(0) is 79% from < k < 1/a and 21% from k > 1/a. Note, 
however, that for the exponentially decaying BEC (7), ip(x) oc (/c|a;|) 1/ i i f I/ (/c|a;|) — > 2 u ~ l T(v) 
cis Ki I x I y for v > 0, which can become arbitrarily large owing to the simple pole in T(v) 
at v = 0. However, for v = 0, f/'f^) oc ^(^l^l) ~~ ► — ln(/c|x|) 0. This logarithmic 

divergence at x = is a direct consequence of the large values of the momentum k, which 
results in the logarithmic divergence of the integral (7). This differs somewhat from the 
study [30| that suggests that, in general, it is the contribution of waves with very small k 
that is important for the accurate determination of i/j(x) in the center of the localized BEC. 
It may be, however, that the momentum k that appears in the Fourier integral (Eq. (6) in 
the first of Ref. 30) may not be so directly connected with the variable k that appears in 
the BEC density (Eq. (8) in the first of Ref. (30)) via the stationary, long-time momentum 
distribution T>(k) and the localized function 4>k{z) of the plane-wave component e lkz . 

It is interesting that the behavior of the BEC density at large distances is determined by 
momenta near the high-momentum cutoff k c owing to the localization of the independent k 



waves 



30l | . For £j n > or, the Lyapunov exponent has a finite lower bound that leads to a 



BEC density that is exponentially localized. On the other hand, for < an, there is no 
such finite lower bound and so the localization is algebraic 30|. It is important to remark 
that the Bogoliubov inequality requires a limit point (or accumulation point) of plane wave 
momenta at k = for the existence of a BEC for D < 2 and that the high-momenta waves 
do contribute to the limit point and so to the removal of the l//c 2 -singularity. 



VI. THE GROSS-PITAEVSKII EQUATION 

The numerical results presented in Ref. (30) for the dynamic behavior of the BEC are 
based on the time-dependent Gross-Pitaevskii equation (GPE) [341 ] . The GPE represents 
a mean-field description of the ground state and it is obtained by finding an extremum (a 
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minimum) of the energy as a functional of the BEC wave function 35]. The time- independent 
GPE for a conserved number of particles corresponding to the Hamiltonian (1) is 

-h 2 r 

— VVW + V ext (r)t/j(r) + V>(r) J dr'[V(r - r') + V(r' - r)] |^(r') | 2 - //^(r) = 0, (45) 

where \x is the chemical potential. 

It is important to remark that the dynamical symmetry-breaking term (5) that requires 
macroscopic occupation of many single-particle momentum states is determined solely by 
the two-body interaction potential V(r — r') and not at all by the external potential V ex t(r). 
Accordingly, the consistency proviso that makes the symmetry-breaking term (5) vanish, 
viz., that the operator y^(r) be orthogonal to both ip(r) and x( r ); does not follow from the 
GPE (45). In fact, for the ground state, our consistency proviso, which requires the last 
two terms in (45) be orthogonal to the operator <^(r) , requires, therefore, that the sum 
: Y^-'V 2 ^(r) + V ex t(r)ift(r) in the GPE (45) be also orthogonal to ${v). The kinetic energy 
term is certainly orthogonal to v^(r) since the kinetic energy is diagonal in the single- 
particle momentum representation; however, the term V ext (r)ip(r) in (45) is not diagonal in 
the single-particle momentum representation and, therefore, the Fourier components V ext (k) 
of V ext (r) must be in the set {k'} of condensate vectors, that is, k G {k'}. 

For instance, for the BEC in the one-dimensional harmonic trap of Sec. II B, 

V ext {x)^{x) = J^Y, eiqX Y. V ™tm q -k, (46) 

V L Q k 

where L is the length of the one-dimensional "box." Our consistency proviso requires that 
(46) be orthogonal to y3^(r); therefore, not only q G {k'} but also k G {k'} since otherwise 
would vanish since ^ q -k is nonzero only for (q — p) G {k'}. Accordingly, the Fourier 
coefficient V ext (k) of V ext (x) cannot have any nonzero Fourier components outside of the 
single-particle momenta that constitutes the condensate, viz., {k'}. The Fourier expansion 
of the harmonic trap is given by 

2 71-2 , a v^/ ^ n cos(k nx) 
x =^2+4E(- 1 ) p-a -n/k <x<n/k (47) 

with nonzero coefficients for k = nko, n = 0, ±1, ±2, • • •, which are the same single-particle 
momenta with macroscopic occupation of the condensate wave function ipo(x) given by (15). 

It should be noted that we are considering the restrictions placed on BECs by the Bogoli- 
ubov inequality at finite temperatures. The determination of the condensate wave function 
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requires us to minimize the Helmholtz free energy with respect to the condensate wave 
function for fixed density and temperature. However, such thermodynamic potential is not 
available. Therefore, the study of BECs given by the GPE does not suffice since the GPE 
describes the ground state (zero temperature) of bosonic systems when all the particles are 
in the condensate. Nonetheless, we have shown that the symmetry breaking term (5), to- 
gether with the GPE, imposes conditions on the momenta of the Fourier coefficients of the 
external potential. 



VII. SUPERSOLID BEC 



The analysis of the possible existence of BEC for D < 2 in the previous sections was based 
on two-particle, local interactions. It was shown that for periodic potentials, the condensate 
wavefunction can be of the Bloch form only for D = 3. Clearly, local interparticle potentials 
cannot give rise to a BEC of the Bloch form for D = 2 since the l//c 2 -singularity in the 
Bogoliubov inequality cannot be removed and still preserve the Bloch form for the BEC. 

The k 2 behavior of the double commutator (21) follows from the kinetic energy term of 
the Hamiltonian since local interparticle potentials V do not contribute to the Bogoliubov 
commutator, viz., {[\C, V], C^}) = 0. It is interesting that the latter is not the case for 
nonlocal potentials |36[. If, for instance, the two-particle potential is a sum of a local 
and a nonlocal potential, then the former potential does not contribute to the Bogoliubov 
commutator while the latter does and if the decay of the nonlocal potential with distance 
is sufficiently slow, then ([[C, V], C* 1 "]) oc k 2 ~ e with e > as k -»■ Q- Therefore, the 
symmetry breaking term (5) allows a BEC of the Bloch form for D = 2 in the presence of 
an infinitely long-range nonlocal potential between the condensate atoms. 

Recently, a Dicke quantum phase transition was realized in an open system formed by 
a BEC coupled to an optical cavity that gives rise to a self-organized supersolid phase 



37j |. It is interesting that the phase transition is driven by infinitely long-range interactions 



between the condensed atoms. The analogy of that work to the Dicke model is based on the 
interaction Hamiltonian that gives rise to a coupling of the pump and cavity fields to the 
zero-momentum states of the atoms to the symmetric superposition of atomic states that 
carry an additional unit of photon momentum. This is quite analogous to our dynamically 
generated symmetry breaking term that allows condensation in atomic states that are integer 
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multiples of a given condensate momentum. 



VIII. SUMMARY AND CONCLUSION 

In random potentials, the underlying mechanism for AL is the suppression of particle 
transport due to destructive interference. The intriguing question is if such mechanism is 
undermined by the presence of interparticle interactions. For bosons, we have seen that the 
generation of periodic BECs is a direct consequence of the dynamical symmetry-breaking 
term in the Hamiltonian that results from the macroscopic occupation of just two single- 
particle momentum states, viz., k = 0, k . However, such types of BECs, albeit allowed 
for D = 3, violate the Bogoliubov inequality for D < 2. Therefore, the presence of atom- 
atom interactions requires that ko — > and thus k = becomes an accumulation point 
of condensates. Note that for D — 1, that suffices to remove the l//c 2 -singularity in the 
Bogoliubov inequality. However, for D = 2, the removal of the l//c 2 -singularity requires 
augmenting the set {k } of condensed states so that the removal of the 1 //^-singularity 
occurs for all approaches to the origin k = 0. 

Finally, the existence and nature of BECs for D < 2 for systems with noninteracting 
atoms in a disordered potentials or for systems in a "box" with only atom-atom interactions 
is precisely the same. Therefore, it may be difficult to discern the individual contribu- 
tions to localization owing to the strength of the disordered potential or the strength of the 
interparticle potentials. Experiments with ultracold atomic gases, where the effect of the dif- 
fering interactions may be easily controlled, will certainly help resolve these very important 
theoretical questions. 
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